Decoherence of Fock States Leads to a Maximally Quantum State 
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We consider the Wigner function evolution of Fock states \n) linearly coupled to a Markovian 
bath of oscillators. In the absence of environmental coupling, apparent "quantumness" increases 
with n, but the presence of any environmental interaction causes high-n states to lose their quantum 
features more rapidly than low-n states. Using the negative volume of the Wigner function as a 
metric [l|, we observe a time-dependent quantumness peak across the eigenstates. 

PACS numbers: 03.65.Yz 
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The macroscopic world is constructed of fundamentally 
quantum objects, but the observation of quantum effects 
requires either very small or very cold systems. That is, 
quantum effects are most easily visible in systems suffi- 
ciently isolated from their environments to prevent de- 
coherence. It is generally more difficult to adequately 
protect larger systems from decoherence; this difficulty 
gives rise to the classical behavior of the macroscopic 
scale [2!] . Contemporary investigations of the transition 
from quantum to classical are motivated by fundamen- 
tal considerations as well as practical issues of control 
and engineering. The transition is being investigated in 
a wide array of systems, including nanomechanical sys- 
tems J3|, y , mesoscopic systems such as Josephson junc- 
tion devices [5[, as well as cavity-QED systems [6(. 

This quantum-classical transition is understood to be a 
multi-parameter transition, controlled by internal system 
parameters and behaviors as well as (a) the relative size 
of h compared to the characteristic action of the system, 
(b) the strength of the system-environment interaction, 
and (c) the temperature and other characteristics of the 
environment, leading to a rich landscape of possible be- 
haviors. 

In this paper we argue that the interplay between the 
relative effects of decoherence and increasing quantum 
number leads to nonmonotonicity in the "quantumness" 
of a unipartite system. We illustrate this in the behavior 
of Fock states \n) for a harmonic oscillator coupled to a 
environment modeled as a Markovian bath of oscillators 
and demonstrate the existence of a time-dependent quan- 
tumness peak. Signatures of this effect have been seen 
in recent experimental work mapping the decoherence of 
Fock states [T], |8| , although the significance had not been 
previously realized. 

In the interaction picture, the time evolution of the 
density matrix p for a quantum harmonic oscillator cou- 
pled to a continuum of oscillators in the Markovian ap- 
proximation and with appropriate assumptions about the 
spectrum of the environment is governed by the master 
equation [9[ 



Here the dot represents the time derivative, the Lindblad 
superoperator is defined as L[0]p = 20pO^ — O^Op — 
pO^O, a^ and a are the raising and lowering operators 
for the harmonic oscillator, respectively, 7 represents the 
degree of coupling of the oscillator to its environment, 
and N corresponds to the mean number of thermal pho- 
tons in the bath. 

In our analysis below, we consider the Wigner function 
W{a), a phase space representation of the state of the 
system, which is defined most easily (if circuitously) via 
the Fourier transform of x(-^)i the Wigner characteristic 
function [9|: 



W{a) 



— x(A)exp(A*a-Aa*), (2) 
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where a = (x + ip)/\2, in terms of the phase space 
quadratures x and p; and \ — {x + iy)l\/2 is defined 
on the Fourier transform of phase space. The func- 
tion x(A) is written in terms of the density matrix as 
x(A) = Tr [pexp(Aa''' — A*a)] . The master equation ([T]) 
can then be rewritten [lO| as a diffusion equation for the 
characteristic function, %: 



7 
^ = -2 






(3) 



Restricting our attention to systems initialized in the nth 
Fock state, the solution to Eq. ([3]) can be written [10| 

X„(A,i)=i„(|Ape-^*)e^l^l^(i-2^('=-^*-i)), (4) 



where L„(x) = Em=o 
mial of order n. 



i-^r 



') is the Laguerre polyno- 



Recall that in the absence of environmental interac- 
tioii, the Wigner function for the nth Fock state is given 
by 



W. 



'„(x,p) = Ljyi^-i^^+P^)L^^2{x'+p^)). 



(5) 
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p^l{NL[a^p+{N + l)L[a]p). 



(1) 



Given Eq. (O, for zero temperature {N = 0), the 
time-dependent solution (j4]) for an environmentally- 
interacting harmonic oscillator becomes a time- 
dependent superposition of non-interacting Wigner 



functions Wm{x,p), where m < n: 
W/„(a;,p,t)=e-"^*f] f ") (e^*-ir-"W™(a;,p). (6) 



m=0 



That is, the Wigner function of a decohering Fock state 
is composed exclusively of diagonal elements. 

To demonstrate Eq. ([6|) , we note that substituting Eq. 
dH) into Eq. ^ yields 



Wnix,p,t)^ J2 



(_g-7t)m /^ 



m=0 



(27r)^?7i! \m 



fm{x,p), (7) 



where 



f^{x,p) = dxdp\X\^"'e^^^^ exp(A*a-Aa*). (8) 
Proposition. We must now establish that 

!,^{x,p) = (2^)2m!^(-l)M "MWfc(x,p). (9) 



fc=0 



Proof. We prove Eq. (|9]) by induction. First, we rear- 
range Eq. ^ to get a recursion relation fn{x,p), yield- 
ing, 



fnix,p) = 

n-l 

-E 

m=0 



2r,I 



(27r)2n! 



Wn{x,p,t) 



(_g-7t)m /^ 



(27r)277i! \m 



frn{x,p) 



for n > and 



/o(a:,p) = (2^)2W^o(x,p,t), 



(10) 



(11) 



which is the base case {n = 0) of Eq. ^ . To prove that 
^ holds true in general, we now show that Eq. © holds 
for the n-f 1th case. We consider fn+i{x,p) by evaluating 
Eq. ([TOl) at t = 0: 



fn+l{x,p) 

n 

-E 



_ (27r)^(n + l)! 
(-1)"+! 

(-1)" /n + 1 



m=0 



(27r)2m! \ m 



Wn+lix,p,t) 



fr,i{x,p) 



(12) 



Because the sum only considers fm{x,p) where m < n, 
we may employ Eq. Q, giving 

/,,+i(i,p) = (2„f(n + ly.l [(-l)"+'»'„+i(i.p,f) 

+ i;i:(-ir"«(":n(";)w.,P)|. (13) 



The term in the square brackets is the n + 1th term of the 
sum in Eq. (J9l); we must demonstrate is that the double 
sum in Eq. (fT3|) constitutes the first n terms. The double 
sum can be rewritten as 



fc=0 ^ ^ 

(-l)"W„+i(x,p). 



(14) 



The second term in Eq. ()14p constitutes the first n terms 
of the sum in Eq. ([SJ evaluated for n-f-1. All that remains 
to be shown to establish Eq. (|9]) is that the first and third 
terms in expression ()14p cancel exactly; that is, we must 
show 

n+l m / I 1\ / \ 

E Ei-^r-rt) nw.(x,;.) = w„,r(x,p). 

r?j.=Ofc=0 ^ / \ / 



The LHS of Eq. (|15p can be rewritten as 

n+l n+l 



Ew™E(-i)' 



k-\-m 



m.=0 fc= 



n + l\ / k 



(15) 



(16) 



but the second sum in Eq. (J16p is the well-known com- 
binatorial identity 



fc= 



kj \a 



(17) 



where 5a. fo is the Kronccker delta [llj; therefore the only 
surviving term in the LHS of Eq. (|16p is yV„+i{x,p), 
thereby establishing Eq. ([9|). ■ 

Inserting Eq. © into Eq. ^ gives 



W„ix,p,t) = X: E(-l)'+"^-™^*Q (7)w.(x,p); 

(18) 
by collecting like Wk{x^p)^ we can rewrite Eq. ffT8|) as 

n 

W„(a;,p,i)= ^ W„(a;,p)( 



771—0 



E/_-| \fc+7rt -fc-yt {n — m). 



(fc — ?n)!(n — A:)! 



(19) 



By changing the index on the second sum to j = k — m, 
Eq. ([TO|) becomes 



W„(a;,p,i) = ^ Wm(a;,p) 



m— 
n— m 



X e-"T* ^ (_l)J(e^*)("-™)-J' f " . "" ) ; (20) 




(a) 



(c) 



FIG. 1; (a) Evolution of the Wigner function representations of a harmonic oscillator initialized in the n = 1,4, and 7 Fock 
states. At 'yt = 0, the strongest quantum features are exhibited by the n — 7 state, but interaction with the environment yields a 
nearly classical distribution by 7i = 0.3. The oscillator initialized in the n — 1 state looks "less quantum" at 'yt = 0, but retains 
its quantum features throughout the evolution, resulting in a comparatively "more quantum" distribution at "/t = 0.3. (The 
fourth quadrant of the xy -plane has been removed in order to emphasize the negative values obtained by the Wigner function.) 
(b) (color online.) Negative volumes of the Wigner function representations for oscillators initialized in the n — 1,4 and 7 Fock 
states as a quantumness metric to quantify the behavior observed qualitatively in (a), (c) (color online.) Comparison of the 
numerically-calculated negative volumes for harmonic oscillators initialized in the 0-20th Fock states for four increasing values 
of jt showing a time-dependent negative- volume peak. The curves from top to bottom correspond to yt — .15, .20, .25, and 
.30, respectively, and the dashed lines indicate the systems initialized in the n = 1,4 and 7 Fock states. 



but the second sum in Eq. ()20p is simply the expansion 
of (e'^* — 1)"-™, and hence Eq. ([6]) is recovered. 



Fig. 1(a) shows plots of the Wigner functions for 



zero-temperature harmonic oscillators initialized in the 
n = 1,4: and 7 Fock states for three increasing values of 
'yt. (Note that allowing the system to evolve for a longer 
time i at a given degree of coupling 7 is equivalent to al- 
lowing the system to evolve for a shorter time at a higher 
degree of coupling — increasing jt captures either of these 
behaviors.) When 7^ = 0, the system initialized in the 
n — 7 Fock state exhibits strongly quantum features: 
its Wigner function has large-amplitude oscillations and 
contains large regions of negative quasi-probability, an 
obvious indicator of non-classical behavior. In these re- 
spects, the system initialized in the n = 4 Fock state is 
comparatively "less quantum" at "ft = 0, and the system 
initialized in the n — 1 Fock state even less so. 

At 'yt — 0.15, the oscillation amplitudes of the sys- 
tem initialized in the n = 7 Fock state are greatly re- 
duced due to environmental interaction. The influence 
of the environment is present in the other two systems 
as well, but the effect is less pronounced. By 'yt = 0.3, 
the systems initialized in the n = 4 and 7 Fock states 



resemble classical orbits. Though it began as the "least 
quantum" state, by 'yt = 0.3, the Wigner function of the 
system initialized in the n — 1 Fock state has retained the 
comparatively "most quantum" oscillations and negative 
quasi-probability. 

A numerical measure of quantumness for unipar- 
tite systems is necessary in this context to compare 
behavior. There are various ways of parameterizing 
nonclassicality — it is typical to define a coherent state 
with minimum uncertainty as the most classical system 
and to consider the quantumness of a state in question to 
be the minimum distance to the coherent state, using a 
distance metric based on the trace, the Hilbert-Schmidt 
distance, or a similar metric 13- It is also possible to 



characterize quantum states based on the properties of 
their Wigner function representations, an attractive op- 
tion in light of recent successes reconstructing the Wigner 
function experimentally |13l . Il4l | . Naively, one might con- 
sider the purity, V, defined in phase space as 



V ^2t: 



dxdp[W{x,p)Y 



(21) 



to be a good indicator of quantumness. However, at zero 
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of decoherence for larger n states is higher. Fig. 2(a) 



n= 1 
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FIG. 2: (a) Peak quantum eigenstate as a function of 7i. (b) 
(color online.) Time-dependent occupation probabilities for 
oscillators initialized in the n = 1 (left) and n = 3 [right) Fock 
states. The solid, dashed, dotted, and dot-dashed lines are 
the occupation probabilities for the n = 0, 1,2 and 3 states, 
respectively. These curves strongly resemble those that have 
been determined experimentally [7|, [a| . 



temperature the time-asymptotic distribution for any ini- 
tial Fock state is |0). The purity of |0) is unity, but its 
Wigner representation >Vo(a;,p) is gaussian, resembling a 
very classical distribution. 

Kenfack and Zyczkowski ll[ suggest that the nonclassi- 
cality, or quantumness of a system can be captured using 
the negative volume of the Wigner function, here denoted 
77 and given by 



ri{n,t) = dxdp 



\Wn{x,p,t)\-Wnix,p,t) 



(22) 



Eq. (|22|) can be numerically integrated to obtain the neg- 
ative volume of the zero-temperature Wigner functions. 



In Fig. 1(b) the negative volumes for the oscillators ini- 
tialized in the n = 1,4 and 7 Fock states are plotted as 
a function of jt. 

Fig. 1(c) shows the negative volume for the 0-20th 
Fock states at four different values of 'yt. For jt = 0, rj 
increases monotonically across the Fock states (approxi- 
mately as ^/n), as reported by Kenfack and Zyczkowski 
m. As the system evolves, rj obtains a finite-n peak 
value, because though oscillators initialized in higher 
Fock states are "more quantum" at the outset, the rate 



shows the peak quantum state as a function of •yt. The 
existence of a quantumness peak can be understood pre- 
cisely as resulting from these two competing trends: (a) 
the tendency for the fringing behavior (and therefore the 
negative volume) to increase with initial n, and (b) the 
robustness against decoherence to decrease with n. This 
behavior is arguably independent both of the measure of 
quantumness and of the system under consideration. 

As mentioned above, the Wigner function for the in- 
teracting oscillator is composed exclusively of diagonal 
elements in the harmonic oscillator basis. Its density 
matrix representation can therefore be written /5„ = 
Y^2=oPk{t)\k){k\, where Pk{t) is the probability of ob- 
serving the system in the fcth eigenstate at time t. Since 
the density matrix is composed of mutually orthogonal 
states, we can use the phase space trace to determine the 
probabilities pk (t) [15| . For an oscillator initialized in the 
nth Fock state 



Pk{t) 



^ e""'^* 



{e'' ~ 1) 



i — k 



k< 



(23) 



Because the density matrix contains no off-diagonal el- 
ements, these probabilities are a complete description of 
the quantum state. The behavior oi pk{t) for oscillators 
initialized in the n = 1 and n = 3 Fock states are shown 



in Fig. 2(b) The occupation probabilities for decohering 
Fock states up to n = 15 have been determined experi- 



mentally T, 8] . The strong resemblance of the curves in 
Fig. 2(b) to experimental curves [l6| suggests that evi- 
dence for a quantumness peak has already been observed. 

In conclusion, the Wigner function representation of 
decohering zero-temperature Fock states can be used to 
compare the quantumness of systems, using the negative 
volume of the Wigner function as a quantumness met- 
ric. Though systems initialized in higher-n Fock states 
start out as the "more quantum" states, they decohere 
more quickly than lower-n states; the competing trends 
of quantum number and decoherence give rise to a time- 
dependent quantumness peak. The general nature of this 
argument suggests that this quantumness peak may be a 
universal feature of unipartite systems. 
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